Measurement of Linear Stark Interference in 199 Hg 
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We present measurements of Stark interference in the 6 1 So — > 6 3 Pi transition in 199 Hg, a process 
whereby a static electric field E mixes magnetic dipole and electric quadrupole couplings into an 
electric dipole transition, leading to i5-linear energy shifts similar to those produced by a permanent 
atomic electric dipole moment (EDM). The measured interference amplitude, asi = (ami + oe2) 
= (5.8 ± 1.5)xl0~ 9 (kV/cm)- 1 , agrees with relativistic, many-body predictions and confirms that 
earlier central-field estimates are a factor of 10 too large. More importantly, this study validates 
the capability of the 199 Hg EDM search apparatus to resolve non-trivial, controlled, and sub-nHz 
Larmor frequency shifts with EDM-like characteristics. 

PACS numbers: 32.60.+i, 32.80.-t,42.50.Gy 



By mixing opposite-parity states, a static electric field 
E can induce an electric dipole (El) amplitude in a mag- 
netic dipole (Ml) transition; one of the resulting Stark 
interference (SI) effects enables precise measurements of 
atomic parity nonconservation (e.g., Refs. [UH]). Analo- 
gous effects occur when an £ , -field mixes Ml and electric 
quadrupole (E2) amplitudes into an El transition. This 
type of SI, first identified [3] and observed [3] with Rb and 
later calculated for Hg [51 [3] , is relevant to searches for 
permanent atomic electric dipole moments (EDMs) [TJ: 
SI is similar to a finite EDM as both produce i?-linear 
ground-state energy shifts that depend on the electronic 
or nuclear spin. Thus SI is a potential source of system- 
atic error in EDM searches and also serves as a useful 
proxy for EDM effects. 

We have measured SI in the 199 Hg — > 6 3 Pi transi- 
tion, the 254 nm El intercombination line used to search 
for the EDM of 199 Hg [HUH]- To our knowledge, this is 
the first observation of SI in a diamagnetic atom, where 
the ground-state polarization is specified by the nuclear 
spin. The measured SI amplitude agrees with relativis- 
tic, many-body predictions [6], and confirms that earlier 
central-field estimates [5J were a factor of 10 too large and 
thus overestimated SI systematic errors in 199 Hg EDM 
searches. In addition, this study was conducted with the 
199 Hg EDM search apparatus, and demonstrates its ca- 
pability to resolve non-trivial, controlled, and sub-nHz 
Larmor frequency shifts with EDM-like signatures. 

199 Hg has nuclear spin 1—1/2 and thus the 6 3 Pi state 
is a hyperfme doublet with F = 1/2, 3/2. SI on the Q 1 Sq 
—> 6 3 Pi transition leads to a small fractional change in 
the El absorptivity a for both the F = 1/2 and 3/2 
hyperfme components [5j: 

(5a/a) 1/2 = -2(Sa/a) 3/2 = a SI (i ■ E)(k x e) • a (1) 

where the SI amplitude asi = baa + a F i is the sum 
of induced Ml and E2 contributions, and e and k are 
unit vectors for the light polarization and propagation 
direction, respectively, a = 2 < I > is the ground state 
nuclear spin polarization. In the present study, a lies 



along an external magnetic field. 

Both Eq. |l]) and a related expression for Si-induced 
ground-state energy shifts are discussed below. We have 
used these shifts, manifested as i?-correlated modulations 
of the ground-state Larmor frequency Q F , to measure 
asi- The measurements span several vector arrange- 
ments with (e • E)(k x e) • a ^ and together give: 

a S i = (5.8 ± 1.5) x 1(T 9 (kV/cm) -1 

Separate measurements with (e • E)(k X e) • a = give: 

(Sa/a) NuU = (0.6 ± 1.8) x 10" 9 (kV/cm) -1 

For both, the quoted error is a quadrature sum of statisti- 
cal and systematic errors; as shown below, the statistical 
errors dominate. The measured 05/ agrees with the Ref. 
[5] prediction, asi = 8.0xl0~ 9 (kV/cm) -1 within 1.5-cr. 
Moreover, (Sa/a)N U u is consistent with and provides a 
useful check on the expected Eq. vector dependence. 

Eq. is obtained by evaluating the product of the 
El amplitude and the Stark induced Ml and E2 am- 
plitudes. The vector e comes from El, E from the 
Stark mixing matrix, and k x e from the optical mag- 
netic field and electric field gradient that drive the Ml 
and E2 amplitudes respectively. The exact vector ex- 
pression in Eq. is derived rigorously in Ref [3] . The 
absorptivity change 5a due to each hyperfme line must 
be equal in magnitude and opposite in sign, as a nu- 
clear spin-dependent effect must vanish if the hyperfme 
structure is not resolved. Since the absorptivity a of the 
F = 3/2 line is twice that of the F — 1/2 line, this 
requires (Sa/a) 3 / 2 — ~ 0-/^)(^ a / a )i/2 as m Eq. Uj. 

SI produces an energy shift SU in the groundstate 
magnetic sublevels and a change <5r in the photon absorp- 
tion rate T. We can relate i5r to the absorptivity change 
via ST /T = Sot/ a. T and ST have a Lorentz lineshape 
7 2 /(4Aw 2 + 7 2 ) with Aw the detuning from resonance 
and 7 the resonance FWHM, while SU has a dispersion 
shape and is related to ST by the dispersion relation [TU] : 

5U F =h^^ST F = h^^T F (Sa/a) Fl (2) 
If If 
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FIG. 1: (color online). Simplified diagram of the apparatus. 
OT, outer top; MT, middle top; MB, middle bottom; OB, 
outer bottom. The fc-vectors for the middle and outer beams 
point along x and z, respectively. The iJ-field is along ± y 



where the subscript F = {1/2, 3/2} denotes the contribu- 
tion of each hyperfine component individually. The work 
presented here uses a laser [TT] tuned midway between 
the two hyperfine components. In this case, the shift to 
fiz, in an external magnetic field B is given by 
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where A is the hyperfine splitting and T to t — ^1/2 + ^3/2 
is the experimentally relevant total photon absorption 
rate; we have used r 3 / 2 = ZTi/2 an d the fact that SQl 
is twice the shift of an individual magnetic sublevel. The 
Eq. spin dependence allows 5VLl to be modeled as a 
Zeeman shift induced by a virtual magnetic field along 
k x i. Hence, shifts linear in E arise if e, E, and B are 
aligned such that k x e has a finite projection along B 
and e has a finite projection along E. 

Fig. [l] is a simplified diagram of the apparatus. De- 
tailed descriptions appear elsewhere [51 H] . Four isotopi- 
cally enriched 199 Hg vapor cells are placed in a uniform 
B-field. The middle cells have oppositely directed in- 
fields, giving Si-sensitive Larmor frequency shifts of op- 
posite sign. The outer cells, placed at E = inside the 
high voltage (HV) electrodes, are free of SI effects and 
instead allow cancelation of B-field gradient noise and 
tests for spurious HV-correlated B-field shifts. The cell 
end caps are coated with SnO and act as _B-field plates 
separated by 11 mm. The cells contain 475 torr of CO 
buffer gas and have paraffin wall coatings that enable typ- 
ical transverse spin coherence times T 2 of 100-200 sec. 
The cells sit in an airtight, conductive polyethylene ves- 
sel housed inside three layers of magnetic shielding. To 
minimize leakage currents, the vessel is filled with 1 bar 
of N2 which is flushed continuously. 

The interference amplitude and (5a/a)Nuii were each 
measured for \E ■ B\ = 1 and \E ■ B\ = 0. The former 
(latter) used all four vapor cells (the middle cells: here, 



k x B = in the outer cells, leading to inefficient pump- 
ing). FigureJ^shows the relevant coordinate systems. For 
Fig. [2ja), the Eq. ^ angular dependence is sin <j> cos <\> 
where ef> is the angle between e and E. In this case, e 
was set to either </> = 45° or -45° , leading to Si-signals of 
5^l/2 or -50^/2, respectively. Separately, e was set to <f) 
= where the expected shift is zero. For Fig. IJb), the 
angular dependence is cos 2 4>. Here, measurements were 
taken with = and <j> = 90° , leading to expected shifts 
of SflL and zero, respectively. To test the Eq. ([3| probe 
light intensity dependence, complete groups of data were 
taken for T p between ~ 1/600 s _1 and 1/50 s _1 . 

Single experimental cycles (termed scans) use a pump- 
probe sequence; throughout, intensity stabilized 254 nm 
laser light enters each cell with k normal to the preces- 
sion axis along B. During the 30 sec pump phase, syn- 
chronous optical pumping with circularly polarized light 
tuned to the 199 Hg 1 S {F = 1/2) -> 3 Pi(F = 1/2) tran- 
sition builds up spin polarization S in a frame rotating 
about B. During the probe phase, the light polarization 
is switched to linear and the frequency tuned midway 
between the F — 1/2 and 3/2 hyperfine lines. The pro- 
cessing S modulates the light polarization angle at fl^- 
This angle is measured, for each cell, with a photodiodc 
after a linear polarizer. The spin precession is monitored 
for 100-200 sec, after which the pump-probe cycle re- 
peats. The HV is ramped to a new value during the 
pump phase, typically alternating between ± 10 kV. 

The Larmor frequencies are determined by fitting 
the photodiode signals with exponentially decaying sine 
waves. Linear combinations of the four frequencies, vot, 
v MTi v MBi an d vqb are then constructed. Here, OT is 
the outer top cell, MB is the middle bottom, etc. For 
\E ■ B\ = 1, v c = (v M t - v M b) - \ {v t ~ vob) has the 
highest SI sensitivity since it maximally suppresses mag- 
netic gradient noise (through 2nd order). For \E ■ B\ — 
0, the useful combination is v m = {vmt — vmb ). The SI 
signal, Ai/si, is obtained from the HV-correlated com- 
ponent of v c or v m via 3-point string analysis !) . Data 
runs last ~ 24 hours and comprise several hundred scans. 
The run-averaged statistical error for Av$i is set by the 
weighted error of the mean multiplied by the square root 
of the reduced x 2 where typically, \ 2 = 2. T p and <j> were 
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FIG. 2: Coordinate system for measuring (omi + clei) when 
(a) E ■ B = 1 and (b) E ■ B = 0. 
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FIG. 3: Measured 5a /a for (a) (e • E)(k X e) • a ^ and (b) 
(e- E)(k x e) • a = 0. Numbers next to the points are sequence- 
average values for 1/T P in sec. The left-hand insets give the 
dataset-wide +B, —B, and final values for 8a/a. The right- 
hand insets show the measured Avsi versus T p . The solid 
line in (a) is Eq. pj with the measured final value for asi- 



fixed during runs. 4> was measured with calibrated polar- 
izers, known to within ± 2°, and changed daily between 
±45° (0, 90°) for \E ■ B\ = 1 {\E ■ B\ = 0). 

The ratio of residual circular to linear polarization for 
the probe beams, set mainly by cell wall birefringence 
and defects in the probe beam waveplates, was typically 
< 0.13. At this level, derived values for Sa/a are negli- 
gibly impacted; e.g., a ratio of 0.3 produces a 2% error. 
The finite circular polarization, however, generates HV- 
independent vector light shifts (oc k ■ B) whose fluctua- 
tions, due to scan-to-scan changes in k ■ B or e, can lead 
to excess noise. For moderate probe intensities and the 
nominal 90° ±0.5° angular alignment between k and B, 
this noise was often negligible. Achieving similar perfor- 
mance at the highest intensities, however, required set- 
ting k to within 66 < 0.1° of k ■ B = 0. To this end, for 
each cell, and prior to each high intensity run, k was set 
such that differential shifts due to flipping the probe po- 
larization from right to left circular indicated 86 < 0.1°. 

The probe intensity noise at Ql is typically 1.5 x the 
shot-noise limit [H] . High intensity | E ■ B | =0 runs used 
the normalized difference between the polarizer outputs 



for each cell. This step reduced run errors, on average, 
by 1.6x, or roughly the expected factor of 1.5\/2 ~ 2. 

The measurement used four vapor cells, four elec- 
trodes, two vessels, multiple vapor cell and electrode ori- 
entations, and several permutations of the photodiode 
acquisition channels. We did not find statistically signif- 
icant correlations between Ai/si and these changes. The 
components were altered between groups of 10-20 runs 
termed sequences; between sequences, the paraffin in- 
side each cell was melted and the outer surfaces cleaned. 
Flips involving the vapor cells, electrodes, and vessels 
used nominally identical components. Each sequence in- 
cluded an equal number of Si-sensitive dipole HV runs (+ 
- + - HV sequence) for the two main S-field directions. 
Within a sequence, the HV ramp rate was permuted on 
adjacent runs between (4/n) kV/s where n = 1,2,4,6. 

The measured SI amplitude was concealed by adding 
an unknown, SI- mimicking offset to vmt and vmb |3 
9 . This fixed blind was revealed only after the data 
collection, data cuts, and error analysis were complete. 

The dataset comprised 181 runs: 47 used \E ■ B\ = 1 
while 134 used \E ■ B\ = 0. In each case, a roughly equal 
number of runs involved non-null and null arrangements 
of v = (e • E)(k x e) • er, with 99 (82) total runs taken for 
the former (latter). The statistical error for the entire 
dataset is 0.13 nHz, or within 2x of Refs. [SJ H]. Fig- 
ure [3] shows the resulting sequence-level values for Sa/a. 
For the expected null data in Fig. |3jb), v ~ 0. Hence, 
Eq. (§ gives Sa/a ~ independent of Az/gj. To avoid 
this artificial zeroing, Sa/a in Fig. [3jb) was calculated 
with v set to the maximum allowed by \E ■ B\; hence, the 
Fig. j3^b) central values are upper limits. 

The sequence values are divided into one point for each 
B-field direction; each point is the weighted average of 
the relevant runs within the sequence. Numbers next 
to the points are sequence-average values for 1/T P . In 
all cases, the +B and — B data are in good agreement. 
The weighted average of all the +B and — B data for the 
non-null and separately, null vector arrangements also 
agree within l-cr. For both Fig. [3]ja) and[3jb), Sa/a 
is also constant (within errors) over the ~ 10 x change 
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FIG. 4: 5a/ a versus <f> when E-B = 0. The inset shows similar 
measurements for E • B = 1. Solid lines give 5a/ a versus <j> 
predicted by Eq. with the measured (a M i + asi)- 
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TABLE I: Systematic 


error budget (10 


10 cm/kV). 


Source 


E-B = 


E ■ B = 1 


Leakage Currents 


3.29 


9.68 


Parameter Correlations 


9 Q1 

z.yi 


1D.D 


Charging Currents 


0.92 


0.85 


E 2 Effects 


0.62 


1.32 


Vector Alignment 


0.13 


0.13 


Quadrature Sum 


4.53 


19.3 



in r p and thus consistent with Eq. (jlj). In contrast, the 
right-hand insets show the expected Avsi dependence 
on r p : the Fig. [3^ a) slope from a least-squares linear fit 
constrained (not constrained) to pass through the origin 
is 63 ± 20 nHz/s (65 ± 14 nHz/s), while Eq. ^ and 
the measured asi give 59 ± 15 nHz/s. For Fig. |3jb), the 
slope for a linear fit constrained (not constrained) to pass 
through the origin is unresolved and equal to 8.5 ± 12 
nHz/s (16 ± 26 nHz/s). For both Fig. ga) andgb), the 
fit intercept is unresolved at the 0.5-ct level. Note that 
when B is flipped, systematics that change sign relative 
to the SI signal can appear in +B 7 —B differences, but 
will cancel in +B, —B averages. Although the data are 
apparently free of such problems, sequence-level values 
were determined from straight +B, —B averages. 

Figure [4] shows Sa/a versus for \E ■ B\ =0. The 
central values for </> = and (j> = 90° are separated by 
> 2-ct and agree with Eq. ([!]). The inset gives similar 
measurements for \E ■ B\ = 1. Although differences be- 
tween cf> — 45°, -45° and are not resolved, da/ a flips 
sign and passes through zero as predicted by Eq. ([!]). 

Using the weighted mean of the sequence values: asi — 
(5.8 ± 1.4 sta4 )xl(T 9 (kV/cm)- 1 and {8a/a) NuU = (0.6 
± l-7 s tat) x 10 -9 (kV/cm) -1 . These final values include 
a numerically modeled, ~ 3% correction due to the op- 
tical rotation-induced, time-varying projection of e onto 
E. Note if \E ■ B\ = 1 and \E ■ B\ = are considered 
separately, the relevant central values agree within 0.6-er. 

Table [I] summarizes the systematic errors. The leakage 
current error was conservatively estimated from worst- 
case scenarios for the single-cell currents: helixes for \E ■ 
B\ = 1 and lines normal to B for \E ■ B\ = 0. For the 
former (latter), the average single-cell current was 0.42 
pA (0.5 pA). In both geometries, the Avsi versus leakage 
current correlation slope was statistically unresolved. For 
\E-B\ = 1, we use a 1/2 turn loop set by the cell geometry 
[5], an effective current of v^xO.42 pA = 0.59 pA (since 
fields in adjacent cells can add or subtract), and divide by 
2 to account for averaging over the uncorrelated current 
paths in the four cells. A similar calculation was used 
for \E ■ B\ = 0, but with the helical loop replaced by a 
V2x0.5 pA = 0.71 pA line current normal to B. 

The parameter correlation error is the quadrature sum 
of 5pi where each Spi is the product of the HV correlation 



for a given experimental parameter and the correlation 
of that same parameter with Ai>si- Specific parameters 
are: the vapor cell spin amplitudes, lifetimes, relative 
phases, and UV transmission; the laser power, frequency, 
drive current, and control voltages; the three axis ambi- 
ent magnetic field; and the B-ficld coil currents (main coil 
and up to three gradient coils). No statistically signifi- 
cant correlations were found. Under these conditions, the 
parameter correlation error is set largely by uncertainties 
in the correlations (and thus affected by the number of 
runs), leading to the larger value for \E ■ B\ = 1. 

The vector alignment error accounts for angular mis- 
alignment between e, k, E, and B and scales as the prod- 
uct of these errors and the measured asi- The E 2 error 
for \E ■ B\ — is a \-a upper limit for {5a/ q)muII multi- 
plied by the measured < 2% iJ-flip asymmetry [9 . The 
remaining Table I entries are given by Eq. (J3J) via errors 
in Ai/g/; the conversion uses an inflated misalignment an- 
gle of <fi — ± 45° and dataset-wide means for T p : 1/130 
sec for \E-B\ = 1 and 1/60 sec for \E-B\= 0. The Av SI 
error analysis techniques are detailed in Ref. [5]. 

Combining the relevant statistical and systematic er- 
rors in quadrature then gives asi and (Sa/a)Nuii for 
\E ■ B\ = 1 and \E ■ B\ = 0. Taking the error- weighted 
mean of these intermediate values: 

a S i = (5.8 ± 1.5) x 10~ 9 (kV/cm) -1 
(Sa/a) NuH = (0.6 ± 1.8) x lO" 9 (kV/cm)" 1 
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